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One of the main quests in quantum metrology, and quantum parameter estimation in general, is 
to find out the highest achievable precision with given resources and design schemes that attain that 
precision. In this article we present a general framework for quantum parameter estimation which 
relates the ultimate precision limit directly to the geometrical properties of underlying dynamics. 
With this framework we present systematical methods for computing the ultimate precision limit and 
optimal probe states. We further demonstrate the power of the framework by deriving a sufficient 
condition on when ancillary systems are not useful for improving the precision limit. 


I. INTRODUCTION 


An important task in science and technology is to 
find the highest achievable precision in measuring and 
estimating parameters of interest with given resources, 
and design schemes to reach it. Quantum metrology, 
which exploits quantum mechanical effects to achieve 
high precision, has gained increasing attention in re¬ 
cent years |ll-[l9l|. where a typical situation is to esti¬ 
mate the value of a continuous parameter x encoded 
in some quantum state px of the system. To estimate 
the value, one needs to first perform measurements on 
the system, which, in the general form, are described by 
Positive Operator Valued Measurements(POVM), {Ey}, 
which provides a distribution for the measurement results 
p{y\x) = Tr{EyPx). According to the Cramer-Rao bound 
in statistical theory [20l - l^ , the standard deviation for 
any unbiased estimator of x, based on the measurement 
results y, is bounded below by the Fisher information: 
5x > , ^ where Sx is the standard deviation of the 

estimation of x, and I{x) is the Fisher information of the 
measurement results, I{x) = v(v\x )( )^[^. 

The Fisher information can be further optimized over all 
POVMs, which gives 


Sx > 


1 


1 

7 ^’ 


( 1 ) 


where the optimized value J{px) is called quantum Fisher 
information [H, [H, [2^ [2^. If the above process is re¬ 
peated n times, then the standard deviation of the esti- 

nJ(px) 

To achieve the best precision, we can further optimize 
the encoding procedures x ^ Px so that J{px) is maxi¬ 
mized. Typically the encoding is achieved by preparing 


mator is bounded by Sx > 




the probe in some initial state po, then let it evolve un¬ 
der a dynamics which contains the interested parameter, 

Po Px- Usually 4>x is determined by a given physi¬ 
cal dynamics which is then fixed, while the initial state 
is up to our choice and can be optimized. A pivotal 
task in quantum metrology is to find out the optimal ini¬ 
tial state Po and the corresponding maximum quantum 
Fisher information under any given evolution (j)x- When 
(px is unitary the GHZ-type of states are known to be 
optimal which leads to the Heisenberg limit. However 
when (px is noisy, such states are in general no longer 
optimal. Finding the optimal probe states and the corre¬ 
sponding highest precision limit under general dynamics 
has been the main quest of the held. Recently using the 
purihcation approach much progress has been made on 
developing systematical methods of calculating the high¬ 
est precision limit [rt-flfll. [l3 | , however how to actu¬ 

ally achieve the highest precision limit is still largely un¬ 
known, as these methods do not provide ways to obtain 
the optimal probe states. Another restriction of these 
methods [3,3 is that they usually restrict to smooth rep¬ 
resentations of the Kraus operators, which is not intrinsic 
to the dynamics. 

In this article, we develop a general framework for 
quantum parameter estimation which relates the ulti¬ 
mate precision limit directly to the geometrical properties 
of underlying dynamics, this provides systematical meth¬ 
ods for computing the ultimate precision limit and opti¬ 
mal probe states without additional assumptions. This 
framework also provides analytical formulas for the pre¬ 
cision limit with arbitrary pure probe states which spares 
the needs of optimization over equivalent Kraus opera¬ 
tors required in previous studies [3, Q. We further demon¬ 
strate the power of the framework by deriving sufficient 
conditions on when ancillary systems are not useful for 
improving the precision limit. 
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II. ULTIMATE PRECISION LIMIT 

The precision limit of measuring x from a set of quan¬ 
tum states Px is determined by the distinguishability be¬ 
tween Pj; and its neighboring states Px+dx [Hi 113 ■ This 
is best seen if we expand the Bures distance between the 
neighboring states Px and Px+dx up to the second order 
of 


^BuresiPxi Px+dx) — ^J{Px)dx , (2) 

where dBuresjpi, P 2 ) = v^2 - 2Fb{pi,P2), here 

Fb{piiP 2 ) = Tr\J pIp 2 pI is the fidelity between 
two states. Thus maximizing the quantum Fisher infor¬ 
mation is equivalent as maximizing the Bures distance, 
which is equivalent as minimizing the fidelity between px 
and Px+dx- If the evolution is given by (jix, Px = 4>x{p) 
and Px+dx = 'Px+dxip), the problem is then equivalent 
to finding out FB[4>x{p),(l>x+dxip)] and the optimal 
p that achieves the minimum. We now develop tools 
to solve this problem for both unitary and non-unitary 
dynamics. 

Given two general evolution (pi and (p 2 of the same 
dimension, we define the Bures angle between them as 
B{pi,p 2 ) = maxpCOS“i[FB( 0 i(p),</> 2 (p))], ^is gener¬ 
alizes the Bures angle on quantum states [^. From 
the definition of the Bures distance it is easy to see 
Uli+Xp dBy^res\^x(^p) T 4^x+dxi_p)\ — 2 2- COS B (^(px j px+dx) ^ 
thus from Eq.Q we have 

Ifi / \1 !• ®[I COS B(^px , Px + dx)] fr,\ 

max J[px[p)\ = hm -— -. (3) 

P dx-^o ax 


The ultimate precision limit under the evolution px is 
thus determined by the Bures angle between px and the 
neighboring channels 


5x > 


\/ 8[l-COS ,<j>x-idx)] r- 

limdx^o 


(4) 


where n is the number of times that the procedure is 
repeated. If px is continuous with respect to x, then 
when dx -+ 0 , B{px, px+dx) B{px,px) = 0 , in this 
case 


jr, / Ni 1. 8[1 COS B(^px I px+dx)] 

max J[px{p)] = hm ^ 

16sin^ 

= 1 ™ -- 

dx^O dx 


= lim 

dx—^O 


AB {pxj Px+dx) 
dx"^ 


(5) 


the ultimate precision limit is then given by 

1 


6x > 


lim. 


dx—^0 ■ 




\/n 


( 6 ) 


The problem is thus reduced to determine the Bures 
angle between quantum channels. We will first show how 
to compute the Bures angle between unitary channels, 
then generalize to noisy quantum channels. 


A. Ultimate precision limit for unitary channels 


Given two unitaries Ui and U 2 of the same dimension, 
since Fb{UipUI,U 2 pU\) = Fb{p,U\U 2 pU\Ui), we have 
B{Ui,U 2 ) = B{I,uIU 2 ), i.e., the Bures angle between 
two unitaries can be reduced to the Bures angle between 
Identity and a unitary. For amxm unitary matrix U, let 
be the eigenvalues of U, where 9j G (—tt, tt] for 1 < 
j < m, which we will call the eigen-angles of U. If ^max = 
9i > 02 >■■■> dm = dmin are arranged in decreasing 
order, then B{I,U) = .^^ben ^max — ^min < 

7r [3ll - l^ . [39l - l4ll| . specifically if 17 = then B{I, U) = 

(A„ax-An.in)t “ ^xnin)t < tt, where Amax(min) IS 

the maximal (minimal) eigenvalue of FI. This provides 
ways to compute Bures angles on unitary channels. For 
example, suppose the evolution takes the form U(x) = 
^^-ixHf^^N (tensor product of for N times, which 

means the same unitary evolution acts on all N 

probes). Then 

B[U{x), U{x + dx)] = B\I, U\x)U{x + dx)] 


It is easy to see that the difference between the maximal 
eigen-angle and the minimal eigen-angle of (^^-'(Htdx^^N 

is ^max ^min — IA^(Aniax | dx \ t Amin | dx \ t). ThuS 
^Q — iHtdx'^^N^ _ ^max~^min ( Ai A^iax I I Amin I ^^1 

Eq.® then recovers the Heisenberg limit 


5x > 


1 1 
\/u(Amax Aniin)l FI 


B. Ultimate precision limit for noisy quantum 
channels 

For a general quantum channel which maps from 
a mi- to m 2 -dimensional Hilbert space, the evolution 
can be represented by a Kraus operation K{p^) = 
Fjp^Fj, here the Kraus operators Fj,l < j < d, 

are of the size m 2 x rn-i, ~ ^mi- The channel 

can be equivalently represented as 

K{p^) = TrE{UEs{\QE){'iE\® P^)uIjs), ( 8 ) 

where jO^;) denotes some standard state of the environ¬ 
ment, and Ues is a unitary operator acting on both sys¬ 
tem and environment, which we will call as the unitary 
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extension of K. A general Ues can be written as 


Ues 


{We 0 -fma) 


'Fi * * • • • *' 

P 2 * * • • • * 

Fd * * 

0 * * 

0 



u 


(9) 


here only the first mi columns of U are fixed and We & 
U (p) {p X p unitaries) only acts on the environment and 
can be chosen arbitrarily, here p> d as p — d zero Kraus 
operators can be added. 

Given a channel an ancillary system can be used to 
improve the precision limit, this can be described as the 
extended channel 


{K ® Ia){p^^) = Y.^Fj ® Ia)p^^{Fj ® Ia)\ 

3 

here p^^ represents a state of the original and ancillary 
systems. Without loss of generality, the ancillary sys¬ 
tem can be assumed to have the same dimension as the 
original system. 

Given two quantum channels Ki and K 2 of the same 
dimension, let Uesi and UeS 2 as unitary extensions of 
Ki and K 2 respectively, we havefs^ 

B{Ki® Ia-,K 2 ® Ia) = min B{Uesi,UeS 2 ) 

UeS1,UbS2 

= mm B{Uesi,UeS2) (10) 

Uesi '■ ’ 

= min B(\Jesi^ UeS 2 )- 

UeS2 


This extends Uhlmann’s purication theorem on mixed 
states [dsjl to noisy quantum channels. Furthermore 
we show in the appendix that B{Ki ® Ia,K 2 G Ia) 
can be explicitly computed from the Kraus opera¬ 
tors of Ki and Ar2[3i|: if Ki{p^) = Yj’l=iFijP^ fI-, 
K 2 {p^) = Yfj=iF 2 jP^Flj, then cosB{Ki G Ia,K 2 ® 
^a') max||p ^||<2 2 ^min(Ariv F here ^Ynin{Fi^ -t- 

ATjy) denotes the minimum eigenvalue of Kw + 
where Kw = Sy > with Wij as the ij-th en¬ 

try of a d X d matrix W which satisfies ||bF|| < 1(|| • || 
denotes the operator norm which equals to the maxi¬ 
mum singular value). If we substitute Ki = and 
K 2 = K^+dx, where K^{p^) = Fj{x)p^F^(x) and 

Kx+dx{p^) = Fj{x + dx)p^Fj{x + dx) with x being 

the interested parameter, then 


cos B{Kx G lA,Kx+dx 0 I a) 


where Kw = '^ijF}{x)Fj{x -\- dx). By substituting 
(t)x = Kx® I A and (t>x+dx = K^+dx ® I A in Eq.®, we 
then get the maximal quantum Fisher information for 
the extended channel K^ 0 I a , 


max J = 


lim max|n^||< 2 ^ T Af^)] ( 12 ) 

dx^c da;2 


In previous studies the operator We in Eq.®, which 
can be arbitrary chosen, was assumed to depend on x 
smoothly [3, Q- As a result, the W in Eo. lfl^ was re¬ 
stricted to unitary operators that depends smoothly on 
X as explained in detail in the appendix IA 21 This restric¬ 
tion was introduced out of computational convenience in 
previous studies, which is not intrinsic to the dynamics. 
The formula here does not have such assumption and can 
be applied more broadly, for example it can be applied 
to the discrimination of quantum channels which is dis¬ 
crete in naturefsj. Also since any W that is not optimal 
gives a lower bound on the precision limit, the formula 
here also provides more room for obtaining useful lower 
bounds. 

The maximization in Eq. (1121) can be further formulated 
as a semi-definite programming and solved efficiently: 
maX|p^||<2 2'^min(Arw T Fw) = 


1 


maximize 

2 ' 



1 hO, 

1 w 

/ } 


Kw + Afjy — tl F 0. 


(13) 


Another advantage of this formulation is that the dual 
form of this semi-definite programming provides a sys¬ 
tematical way for obtaining the optimal probe states, 
which we will show in the next setion. 


III. OPTIMAL PROBE STATES 

Developing systematical methods to obtain the op¬ 
timal probe states are essential for achieving the pre¬ 
cision limit. So far there are only a few cases for 
which optimal probe states are known, mostly for phase 
estimations [llL[l2 .. A systematical way of obtain¬ 
ing optimal probe states for general quantum dynamics 
is highly desired as it will pave the way for achieving the 
ultimate precision limit. We now show how to obtain the 
optimal probe states that achieve the ultimate precision 
limit for extended channels. 

We first provide an analytical formula for calculating 
quantum Eisher information with any given pure input 
states, for both unextended and extended channels, then 
use it to obtain optimal probe states for extended chan¬ 
nels. 

In the appendix we showed that for both unextended 
and extended channels with pure probe states we have [3^ 

FB[K,{p^),K,+dx{p^)] = ||M(p^)||i, (14) 


( 11 ) 
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Fb[{K^ ® Ia){p^^), {K,+d. ® lA)ip^^)] = ||M(p^)||i, 

(15) 

here M{p^) is a d x d matrix with its ij-entry equals to 
Tr[p^F- {x)Fj{x + dx)], and || ■ ||i represents the trace 
norm which equals to the summation of singular values. 
For the unextended channel this formula works for the 
pure probe state p^ = \ips){ips\, while for the extended 
channel although p^^ = \tpsA){tpSA\ is required to be 
a pure state, p^' = TrA{p^^) can be any mixed state, 
which characterizes the advantage provided by ancillary 
systems. 

The above two formulas provide a straightforward way 
calculating the quantum Fisher information with any 
pure probe states, 


J[kAp^)] 

J[iK,®lA)ip^^)] 


lim 

dx—>^0 


8(l-||M(p^)||i) 

dx'^ 


lim 

dx—^0 


8(l-||M(/)||i) 

dx'^ 


(16) 


In contrast to previous studies [3, d, [l^, optimization 
over equivalent representations of Kraus operator is not 
needed in this formulation. In fact ||M(p'®)||i does 
not depend on any particular representation of the 
Kraus operators: if we use a different representation 
of Kraus operators for and K^+dx-, for example 
= J2r'^irFrix) and Fj{x + dx) = ''^IsT's(a;+da:), 

where Uir and Vjs are entries of some unitary ma¬ 
trices U and V respectively, then Tr[p^F- {x)Fj{x + 
dx)] = ^Kx)Fs{x + dx)], thus M{p^) = 

UM{p^)V'^ which has the same trace norm \\M{p^)\\i = 

\\Mip^)\h. 

The optimal probe states for the extended channel can 
then be obtained by minimizing over input states at both 
sides of Eg. dTSl) . 


minFsHKa: ® Ia){p^^), {K^+dx ® lA)ip^^)] 

pSA 

= min||M(p^)||i. 


This can be computed by a semi-definite programming 
formulation for the trace normjd^ as min^s ||M(p'®)||i = 

minimize 

s.t. 


where P, Q are Hermitian matrices. One can verify that 
this is exactly the dual form of the semi-definite pro¬ 
gramming used in Eg. (11311 . From the output p^, we can 
easily obtain the optimal probe state p^^, which is any 
purification with the reduced states equals to the optimal 
. This gives a systematical way to obtain the optimal 
probe states for the extended channel which we demon¬ 
strate through some examples. 


irr(P) + ^Tr{Q) 

P M^p^ 
M{p^) Q 

P^^0,Tr(p^) = l, 


hO, 


(18) 


Consider phase estimation with spontaneous emis¬ 
sion, K^{pq) = Fi{x)poFI{x) + F 2 {x)pq FI{x) , where 

Fi{x) = ^^U{x),F^{x) = (^[j ^^^)p(x), 

U{x) = exp(—f^x). Suppose a pure input 
state p^^ is prepared for the extended channel 
Kx ® I A, then p^ = TrA{p^^), M{p^) = 


/ Trlp^F}(x)Fi(x + dx)] Trlp^F^ (x)F 2 (x + dx)]'\ 
\Tr[p^F^{x)Fi{x + dx)] Pr[p‘^P|(x)P 2 (x-|-dx)] / ’ 
this case the problem can be solved analytically, in the 
appendix we showed that the optimal p^ is given by 
— 0 \ 

, and the corresponding maximal 


^5 ^ 




quantum Fisher information is max J = 


Since the optimal p^ is mixed, an ancillary systems is 
necessary. The optimal input state in this case is any 
pure state p^^ with the reduced state equal to p^, the 
simplest choice of the optimal input state in this case 


is Yi^^lOO) -I- Y^^Y^^|11), which is not a maximally 

entangled state as previously suspected[T3, 1131 . We can 
also use the method to find the maximal quantum Fisher 
information without using ancilla by imposing the condi¬ 
tion that p^ be pure. In that case the maximal quantum 
Fisher information turns out to be p and the optimal 
input state has the form (|0) -I- exp(f0)|l))/-\/2 for some 
e e R[3|. 


For high dimensional systems, we use the CVX 
package in Matlab[43] to implement the semi-definite 
programming of and obtain the optimal input 

states. For example, consider two qubits with inde¬ 
pendent dephasing noises, which can be represented by 
a Kraus operation with 4-Kraus operators: Fi{x) ® 
Fi (x), Fi (x) 0 1^2 (x), F 2 (x) 0 Fi (x), F 2 (x) 0 Fi (x) with 

Fi(x) = ^^U{x), F 2 (x) = y^a 3 C/(x), here U{x) = 
exp(—f^x). It turns out that min^s ||M(p'^)|| can always 
be attained with a pure p^, ancillary systems are thus not 
necessary. In Fig. [T]we plotted the entanglement of op¬ 
timal states, which is quantified by the entropy of the 
reduced single qubit state p"®, at different p. It can be 
seen that there exists a threshold for p: when p exceeds 
the threshold, the optimal state is the GHZ state, which 
is maximally entangled; when p is below the threshold, 
GHZ state ceases to be optimal, with the decreasing of 
p, the optimal state gradually changes from the maxi¬ 
mally entangled state to separable state. Fig. [5] shows 
the quantum Fisher information with the optimal state 
and the separable input state | -f-f), where |-|-) = . 

It can be seen that the gain of entanglement is only ob¬ 
vious in the region of high p, i.e., low noises. Similar 
behaviour is found for more qubits, i.e., there exists a 
threshold for p, above the threshold the optimal state is 
the GHZ state and with the decreasing of p, the opti¬ 
mal state gradually changes from GHZ state to separa¬ 
ble state, and this threshold increases with the number 
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FIG. 1: Entropy of reduced single qubit state, which is used to 
quantify the entanglement of the optimal state, at different rj. 



n 


FIG. 2: Quantum Fisher information with the optimal input 
state and separable input state | + +) for 2 qubits with inde¬ 
pendent dephasing noises. 



n 


FIG. 3: Optimal probe state for 5 qubits under indepen¬ 
dent dephasing noises. The optimal state has the form 
Itp) = aolV'o) + uilV’i) + a 2 |'!/' 2 ), where lipi) denotes the sum¬ 
mation of all basis states with i zeros or i ones, for example 
l-^o) = lOOOOO) -h llllll). 


of qubits. In Figl3] the optimal state for 5 qubits with 
independent dephasing noises is shown, and in Figl4]the 
quantum Fisher information for the optimal state, GHZ 
state and the separable state are plotted. We also cal¬ 
culated the optimal state for 10 qubits with 5 of them 
under independent spontaneous emission and 5 of them 
as ancillary qubits, in Figl5]we plotted the Fisher infor¬ 
mation for the optimal state, GHZ state and separable 
state. 


FIG. 4: Quantum Fisher information for optimal probe states, 
GHZ state and separable state for 5 qubits under independent 
dephasing noises. 



n 


FIG. 5: Quantum Fisher information for optimal probe states, 
GHZ state and separable state for 10 qubits with 5 of them 
under independent spontaneous emission and 5 of them as 
ancillary system.. 


IV. WHEN ANCILLARY SYSTEMS ARE NOT 
HELPFUL 

The formulas developed here not only provides system¬ 
atical methods to compute the ultimate precision limit 
and optimal probe states, but also have wide implica¬ 
tions, which we will demonstrate by deriving a sufficient 
condition on when ancillary systems are not useful for 
improving the precision limit. 

We have seen that in the spontaneous emission case an¬ 
cillary system helps improving the precision limit while 
in some other cases—known examples include unitary, 
classical, phase estimation with dephasing and lossy 
channels (la 1^ . ancillary systems do not help. For a gen¬ 
eral channel it is usually difficult to tell whether ancillary 
systems can help improving the precision limit or not. 
Previously this problem was usually studied case by case 
by comparing the maximum quantum Fisher information 
for the unextended and extended channels. The formulas 
developed here provide a more direct way: from Eq. dm) 
it is obvious that ancillary systems do not help improv¬ 
ing the precision limit if and only if minpS ||M(/9'^)||i can 
be reached at a pure state , which can be checked 
by using the semi-definite programming to obtain op- 
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timal and see if ^ more easily verifi¬ 

able sufficient condition is as following: given a chan¬ 
nel K^{p) = if FHx)Fj{x + dx), 

^ ^ i, j ^ d can be simultaneously diagonalized, then an¬ 
cillary systems do not help improving the precision limit. 
As if there exist a basis such that F- {x)Fj {x + dx) are 
all diagonal, then in that basis only the diagonal entries 
of p^ enter into M{p^ ){a.s the entries of M{p^) are of 
the form Tr[p^F^(x)Fj{x + dx)] which only depends on 
the diagonal terms of p^ if F-{x)Fj(x + dx) is diago¬ 
nal) , other entries can be chosen freely. This means that 
miups |lM(p'®)||i can be achieved by all p^ with the op¬ 
timal diagonal entries, pf^ = ai, this always includes a 
pure state |^s) = hence ancillary system is 

not necessary to achieve miUpS ||M(p‘®)||i for such chan¬ 
nels. This sufficient condition is satisfied by unitary, clas¬ 
sical and phase estimation with dephasing channel, and 
many other channels that have not been categorized be¬ 
fore, for example phase estimation with noises along the 
X and Y directions satisfies this condition, as it can be 
represented with two following Kraus operators Fi (x) = 

\J exp(—and F 2 {x) = ^^(^2 exp(— 

one can easily check that they satisfy the condition thus 
ancillary system does not help improving the ultimate 
precision limit for this channel. 


V. CONCLUSION 

In conclusion we presented a general framework for 
quantum metrology which provides systematical ways of 
obtaining the ultimate precision limit and optimal input 
states. This framework relates the ultimate precision 
limit directly to the underlying dynamics, which opens 
the possibility of utilizing quantum control methods to al¬ 
ter the underlying dynamics for better precision limit ( 3 ^ . 
The tools developed here, such as the generalized Bures 
angle on quantum channels that can be efficiently com¬ 
puted using semi-definite programming, are expected to 
find wide applications in various fields of quantum infor¬ 
mation science. 


Appendix A: Appendix 

In the appendix we give detailed derivation on the Bu¬ 
res angle between extended channels, the analytical for¬ 
mula to compute quantum Fisher information with arbi¬ 
trary pure input state and two examples with dephasing 
noise and spontaneous emission respectively. 


1. Bures angle between extended channels 


and K 2 ® I A can be computed from the Kraus operators 
as following 

cosB{Ki® Ia,K 2 ® Ia) = max 

||w||<i 2 


i.e., 


minFB[(A:i (g) Ia){p^^), {K2 ® Ia){p^'^)] 

pSA 

1 , (Al) 

here Kw = j the ij-fh entry of d x d 

matrix W and the optimization is over all d x d matrices 
W with operator norm ||IK|| < 1. During the proof the 
analytical formulas that compute quantum Fisher infor¬ 
mation for arbitrary pure input states will be derived. 

As the minimum hdelity can always be achieved with 
pure states, we can assume p^^ = |'0Sa)('0sa|. Denote 
Uesii UeS 2 as the unitary extension of Ki and K 2 re¬ 
spectively, i.e.. 


Kiip^) = TrE{UESl{\0E){0E \® 

K2{p^) = TrE{UES2[\^E){^E\ ® P^)Ues2)^ 


where jO^;) denotes some standard state of the environ¬ 
ment, Uesi and UeS 2 are unitary operators acting on 
both system and environment. The general form of Uesi 
can be written as a pTO 2 “Climensional unitary, here p > d, 
p — d zero Kraus operators can be added. 


Uesi 


(VFi 0 


Fii ^ ^ ^ 

F\2 ^ ^ 

F\d * * ‘ ‘ ‘ * 

0 **■■■* 

0 **■■■* 


Ui 


(A3) 


where Wi € U{p) and only the first uii columns of Ui 
are fixed. Similarly 


UeS2 




F 21 * * * * * * 
F 22 * * * ' ' * 

F2d **•••* 

0 

0 **•••*. 


U 2 


(A4) 


In this section we show that for any two given channels 
which maps from mi- to m 2 -dimensional Hilbert space, 
Ki{p^) = Eti and K2{p^) = 

the Bures angle between the extended channels Ki ® I a 


where W 2 G U{p) and only the first uii columns of U 2 are 
fixed. Note that IFi and W 2 only act on the environment 
and can be chosen arbitrarily. 
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Uia, UeS2 = 


Pt rt 

^11 r ]_2 


* 

* * 




VKrf * * 

* * * 


Pid 0 * 




-^11 ^12 * ■ ■ ■ ^Id ^ * 




(Wl ^ Im,){W2 ^ Im,) 


{Wp®I^^) 


F 21 * * 
F 22 * * 


'F 21 * * 
F 22 * * 

F2d * * 
0 * * 


0 * * 

• *■ 

• * 


F2d * * *‘‘ * 

0 **•••* 


0 +*•••+ 


(A5) 


where Wp = WlW 2 S U(p), Wa is the first dx d subma- and Sj=i with Wij as the jj-th 

trix of Wp, i.e., entry of Wd- Then 


Wp = 


Wd 

* 


* 

* 


(A6) 

I 


min Fb[{Ki (g) /a)(|V'Sa)(V'Svi|), iK 2 0 /a)(|^/’sa)(V'Sa|)] 
lbsA> 

= min FB[TrE{UESl ® ^a(| 0£;)(0£;| (g) \'tjjSA){i’SA\)Ul,g-^ (g) lA),TrE{UES 2 ® /a(|0_e)(0e| g) \i^SA){^SA\)Ul;g 2 ® ^a)] 
IbsA) 

= min maxFBiUESi g> -^a(| 0£;)(0£;| (g) |'!/'Sa)(V'Sa|)C7^si g> 7 a, UeS 2 g> 7 a(| 0£;)(0£;| (g) |'0sA)('0SA|)C7|;g2 g> 7 a), 

|'!/’SA> C^ES2 


the second equality used Uhlmann’s theorem which states 
that 


where \tpi) is an arbitrary purification of pi and the max¬ 
imization runs over all purifications \tjj 2 ) of Thus 

















min Fb[{Ki ® /a)(|'0sa)(V’Sa|), {K 2 O ^a)(|V'Sa)('0Sa|)] 

Wsa) 

= min mayLFslUESl ® ^a(| 0£;)(0£;| (g) |'0 Sa)('0Sa|)C^1;si <g> -^A, t^_ES 2 <g> -fA(|0E)(0£:| (g) |'0 Sa)('0Sa|)C^1;S2 ® 

IV'sa) UeS2 

= min maxFB[|0£;)(0£:| 0 |V'sa)(V'Sa|, U\.g-JJES 2 0 -fA(|0£:)(0£:| 0 |'!/'Sa)(V'Sa|)C^es 2 ^£^‘Si ® /a] 

I'i/’Sa) t^ES2 

= min ma.x\{tj}sA\{ 0 E\Ul;g-^UES 2 ^ lA\0E)\'fpSA)\ 

IV'Sa) [^ES2 


= min max |(i/;saK0£;| 

|'0Sa) 


K 


* * 
* * * 


= min mayiKtpsAlKwd Ia\'4’sa)\ 

\i>SA) Wd 

= mm max | Tr J | 

p® Wd 

= min max |Tr(^ p®f/.F 2 j ) | 

P® VKd ^ 


0 -fA|0_E)|^SA)| 


(A8) 


= min max | ^ Tr(p®F/.F 2 j) | 

VKd 

= mm max | T r [Wj M{p^)]\, 


where p'® = TrA{p^^) and M{p^) is a dx d matrix with 
its ij-entry equals to Tr{p^Fl^F 2 j). As Wd is the first 
dX d submatrix of Wp S U{p), \\Wd\\ < 1, here |||| denotes 


the operator norm which equals to the maximum singular 
value. Conversely any Wd satisfies \\ Wd\ \ < 1 can be 
imbedded as a submatrix of a unitary 4^ , thus 


min ^b[(A:i0/a)(|^/’sa)(iAsa|),(A:2 0 Ja)(|V'Sa)(V'Sa|)] = min max \Tr[WjM{p‘- 

lbsA> pS ||lVd||<l 


r 


(A9) 


Let M{p^) = UDV be the singular value decom¬ 
position of M{p^), where U,V G U{d) and D = 
diag{si,S 2 , - ■ ■ ,Sd), then max|ny^||<i |rr[lVjM(p'5)]| = 
ll-^(P’^)lli = til® equality is achieved with 

Wj = V^W. This gives an analytical formula to com¬ 
pute the fidelity with any given pure input state 


we get the analytical formula for quantum Fisher infor¬ 
mation of the extended channel with any pure input state 

IV'sa), J{Ka; 0 Ia{p^^)) = linida;^o ^^^"''^ 2 ^° where 
P‘5 = TrA{\i^SA){'^SA\)- Following the same line of argu¬ 
ment, it can be shown that without auxiliary system the 
same formula holds 


FB[(Ari 0 /a)(|V'Sa)(V’Sa|), {K 2 0 I a) Hi’S a) {'ip sa\)) 
=\\M{p^)h 


1- 


If we substitute ATi, K 2 with and K^^dx, w® 
get FB{px,Px+dx) = ||M(p5)||i, where p^ = ® 

^a(|'0Sa)(V’Sa|) and p^, -\-dx — ^x-\-dx <8 /a(|V^sa)(^sa|), 
then using the connection between the Bures distance 
and the quantum Fisher information 

^ Bures ip X ; Px-\-dx ) — jJ{px)dx^, (AlO) 


Fb [Kx (I ^/’s) (^/’s I), Kx+dx (I ^s) (^s I)] 
= l|M(p^)||i, 


(All) 


and J{Kx{p^)) = limda^-s-o , just in this case 

P^ = \'ips){'ips\ is not the reduced state, but the pure 
input state \'ips){'i’s\- 

Note that when Wj = rr[WjM(p5)] = 
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||M(/9‘^)||i is a positive real number, thus 

max \Tr\WjM{p^)\\ 

\\Wci\\<l 

= max Re{Tr\WTM(p^)]\ 

= I max^ \Tr{WjM{p^) + [WjM{p^)n 

||Wd||<l z 


(A12) 


optimal point that achieves the optimal value can change. 
If we substitute Ki with Kx and K 2 with K^+dx, then 
we get the formula to compute the Bures angle between 
the extended channels 

cos B(Kx 0 I A, Kx+dx 0 I a) 

= imn Fb[{Ki Ia){p^^),{K 2 <^ Ia){p^'^)] 

(A14) 

which then gives the maximal quantum Fisher informa¬ 
tion for the extended channel Kx 0 I a 


then mini^g^) FB[(iv:i 0 /a)(|V’Sa)(V^sa|), (^12 0 

lA){\i’SA){ipSA\)] = mmpsma.xiiw^ii<i^Tr[p^{Kwa + 
where p^ = TrA{\ilJsA){'4’SA\) and 
with Wy as the ij-th 
entry of Wd- For the extended channel since all the 
reduced states p^ forms a convex set, {IFd|||Wd|| < 1} 
is convex and compact and ^Tr[p^{Kw^ + 
a bilinear function of entries of p^ and IF^, thus from 
Sion’s theorem[^. the sequence of min-max can be 
exchanged, so 

= max min ^Tr[p^{Kw + K\^)] (A13) 

= max “1“ 

|| W ||<1 2 

here \m.in{Kw + k\^) denotes the minimum eigenvalue 
of Kw + K\y. We used a different symbol IF after the 
exchange of min-max to emphasis that although the opti¬ 
mal value is not affected by the exchange of min-max, the 


m.axJ{Kx ®Ia{p^^)) 

pSA 


= lim 8 

dx—^0 


= lim 
dx—>-0 


1 - cos B{Kx 0 I A, Kx+dx 0 I a) 
dx'^ 

8(1 — max||w||<i \^miu{Kw + ^vv) 
dx"^ 


(A15) 


where Kw = {x)Fj{x + dx). 

2. Connection with previous studies 

The above formula includes previous studies as a spe¬ 
cial case. Previous results as in [3, Q state that for an 
extended channel Kx 0 I a the maximal quantum Fisher 
information is given by 


d 

maxJ = 4 min || F|(x)Fj(x)|| (A-16) 

where the minimization is over all smooth representations 
of equivalent Kraus operators of the channel Kx- Note 
that this can be equivalently written as 


J 


max J = 4 min || FUx)Fj{x) 

{F^{x)} ^ 


t=i 

d 


= 4 min || N lim 


(fJ( a; + dx) - Pj(x)) [Fj(x + dx) - Fj (x)) 


dx 


dx 


= 4 min 

{F(a:)} 


, 2d - (-^J i^)Pj + dx) + (a: + dx)Fj (a:)) 


dx'^ 

rt p. I _L pt ! 


= 4 


2 - max^^^(,^)j ^min^j^i{Fj{x)Fj{x + dx) + Fj{x + dx)Fj{x))] 


dx"^ 


(A17) 


where the optimization is over all smooth representations equivalent Kraus operators are represented by Fj (x) = 
of equivalent Kraus operators. In previous studies the 
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sf=i and Fj{x + dx) = + 

dx)Fi{x + dx), where ujji{x) is ji-entry of We{x) S U{d), 
and We{x) is required to be smooth with respect to x. It 
is easy to see that in this case Ea. (IA17|) is just a special 
case of Eq. (IA15I) when W is restricted to taking the form 
wI;{x)We{x + dx). 

We provide an example showing that the optimal W 
in Eg. dAlSI) can be non-unitary. 

Consider two 8-dimensional channels with two Kraus 
operators each. The first channel Ki{dx) has the follow¬ 
ing operators: 


Ell = diag(a, a, a, a, a, a, a, a) (A18) 

Fi 2 = diag(/3, -/3, /3, -/3, ifd, -ij3, ijd, -i/d) (A19) 

The second channel K 2 {dx) has the following: 

F 21 = Ell (A20) 

F 22 = diag(/3, 13, -(3, -13, jd. Id, -jd, -/d) (A21) 


where a = V1 — dx'^ and jd = |fia;| are real. Note that 
these two channels can be regarded as two neighbour¬ 
ing channels parameterized continuously in x, (j}{x) = 
( Ki{x) a; < 0 
1 K 2 {x) a; > 0 ' 


We let W = 
Ally), where ATiy 


, and compute ^XminiKw + 


Wll W12 
W 21 W 22 

= wiiFh F 2 i+wi 2 Fh F22+W21FLF21 + 


W22F12F22. First, note that if W 


1 0 
0 0 ’ 


then 


Now consider general W, since all matrices are 

diagonal in this case, we will only need to look at the 
diagonal elements. The diagonal elements of Kyy are 


a(wiia+wi2/3) +/3(w2ia+W22/3) (A22) 

-b - -b (A23) 

- -b - (A24) 

- - - (A25) 

-b — i -b (A26) 

-b -b* -b (A27) 

- -i - (A28) 

- +i - (A29) 

Thus, Tr{Kw) = Sa^wn, thus Tr{ ^ _ 

8a‘^Re{wii). If Re{wii) < I, the average eigenvalue of 
will be smaller than thus the minimum eigen¬ 
value will also be smaller than a^, which is worse than 


the case of W = 


I 0 
0 0 ’ 


thus not optimal. So for optimal 


W, Wll must be 1. With wn = I, the most general W 


with IIWII < 1 is W = Q 
agonal elements of Kyy are then 


with 0 < r < 1. The di- 


And the diagonal elements of 


± /3^r cos(E) ,a^ ± 0^r sin(0). 


If r > 0, one of them must be less than . 
optimal choice of W in this case is W = g [! , 


So the 
which is 


not unitary. And it is easy to check that the minimum 
value o? for FB[{Ki{dx)®lA){p^^)^ {K 2 {dx)®I a){p^^)] 
in this case can be achieved by preparing the input state 
as the maximally entangled state X]i=i I**)- 
We also note that as the set {IE|||IT|| < 1} is a con¬ 
vex set, this allowed a direct formulation with semi- 
definite programming as stated in the main text. While 
in previous studies W is restricted to be unitary which 
does not form a convex set. To circumvent the diffi¬ 
culty previous study has to resort to the Lie algebra of 
the unitaries and formulated the semi-definite program¬ 
ming there insteadfl^. That, however, comes with a cost 
on the computational complexity, which can be seen by 
comparing the size of the constraining matrices in the 
semi-definite programming: the constraining matrices in 
the semi-definite programming of Eg. (1131) have the to¬ 
tal size of 2d + mi, while previous study needs a size of 
mi + dm9(l^. The difference can be significant when the 
system gets large(note that for generic channels d is in 
the order of 77117712 ). For example for Wqubit system, 
7771 = 7772 = 2^, the difference quickly becomes large 
with the increase of N. 


3. Quantum Fisher information and optimal input 
states with dephasing noise and spontaneous 
emission 


We give detailed calculation for the maximum quan¬ 
tum Fisher information and optimal input states for two 
examples, one with dephasing noise and the other with 
spontaneous emission. 

Consider a channel with dephasing noise 

Kx{p) = U{x){^^^^p+^-Y^a3pa3)U\x), 


where U(x) = exp (—77 
and (73 = 


<X2 = 


0 -i 
i 0 


s [0, 1]. tTl — ^ 

77 e [ 0 , 1 ]. 



In this 


case Fi{x) = \J^^U{x), F 2 {x) = ^-^asU^x). Let 

be a pure input state for the extended channel 0 
I A, and p^ = TrAip^^), then 


M{p^) 

f Tr[p^ fI{x)Fi{x + dx)] Tr[p^ fI{x)F 2 {x + dx)] \ 
\Tr[p^F 2 {x)Fi{x + dx)] Tr[p^F 2 {x)F 2 {x + dx)] ) 




C 
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where C = pf;^exp(—i^) + / 9 I 2 exp(i^), D = 
exp(—z^) — pf 2 exp(i^). Thus 


M{p^)M\p^) 


mil mi2 
m 2 i m 22 


where mu = \C\^ + ^ ^ l-Pp, ^12 = ^^4 [(1 + 

+ (1 - p)C*D], m 2 i = ^^^[(1 + p)C*D + (1 - 
p)CD*], m 22 = |Cp + Denote Si and 

S 2 as singular values of M{p^), then sf+ S 2 = Tr{MM^) 
and S 1 S 2 = |det(M)|, thus 

(si + 52 )^ 


=Tr{MM'') + 2\det{M)\ 


1 + 


1 — ?7 


1 — ?7^ 


+ -^\D\^ + -^\C^ - D^\ (A30) 


2 ' ' 2 ' ' 2 
= (Pn + Pf 2 )^ - 2 pnpf 2 ^^[l - cos{dx)] 


= l-4pnp|2’7^ sin' 


da; 


To minimize Si + 52 ; we just need to maximize pfipf 2 ! 
which is achieved with pfi = P 22 = i- Thus 


COs\^]3(^Kx ^ ^Af d^x+dx ^ 


1 9-2 

1 — ? 7 ^ Sin 


dx 

T’ 


Denote si and S 2 as the singular values of M{p) we can 
similarly get 

(si + S 2 )^ 

=Pii + [v^ + (1 - v)'^]P 22 + 2ppiiP22 cos{dx) (A32) 

+ 2(1 -7?)(|P2i|^ + IPIIP22 + ^226*'*'" - |P2in)- 

First observe that 

|p2l|^ + |pilP22+wLe*""-|p2in 
> ||d2iP +P11P22 + ? 7 pLe*‘^'" - IP2in 
= |piiP22 + wLe*"l 

= P 22 JPu + ‘dpPllP 22 cosdx + p'^P 22 , 


where the equality is achieved with p 2 i =0, so 
(si + 32 )^ 

>Pii + + (1 - vf]P22 + 2?7PiiP 22 cos(da:) 

+ 2(1 - p)p 22 \/p\i + 2 ? 7 PiiP 22 cosdx + p^pl 2 



4 + 2ppiiP22 cos dx + p2p2^ + (1 “ p)p 22 ] 


(A34) 

2 


>1 


I p + 2p^ cos dx + p 2 


(1 + vvr 


+ 


l-d 12 

1 + 0 ? 


we then get the maximal quantum Fisher information 
limrf^^o 8 

sistent with previous studies [l 3 , 113 ■ In this case the 
optimal initial state is any pure state p^"^ such that the 
reduced state p^ satisfies pf^ = P 22 = 5 - The simplest 

choice is the pure state i.e. in this case ancil- 

v 2 

lary system is not necessary, which means that the Bures 
angle for and K^+dx remains the same 


COs[il(Ara;, Kx+dxy\ 

= cos[B{Kx ® I A, Kx+dx 0 I a)] 


= 1 /1 — 34442 


dx 

T’ 


(A31) 


thus the maximal quantum Fisher information under the 
evolution is also p^. This is consistent with previ¬ 
ous studies, however in our framework the optimal state 
comes out from the calculation naturally while in previ¬ 
ous studies it depends on educated guess [l^ 

For a channel with spontaneous emission, the Kraus 
operation takes form Kx{po) = Fi(a;)po+ 4 ^(a:) -f 

+ 2 (x)po+|(a;), where +i(x) = ^ U{x),F 2 {x) = 

0 ^ ) U{x),U{x) = exp(—z^x). Again let p^^ 

be a pure input state for the extended channel 0 I a, 

and p = Trji{p^^), then 


The last inequality is achieved by substituting p 4 i = 1 — 
P 22 into the equation and find the minimum of a single 
variable function, which is saturated with P 44 = 1 +^ 
and P 22 = tT^- T^ie minimum is thus achieved at p = 
0 \ 

I ‘-+v'd ^ ^ which gives the Bures angle 


0 


1 + +5) / 

cos[B{Kx O Ia, Kx+dx ® Ia)] 
p + 2pyyp cosdx + p"^ 1 — p 


(A35) 


( 1 +vvy 


1 + 0 ?’ 


expanding both sides up to the second order of 
dx, we get the maximal quantum Fisher information 

limrf,,^o = (T^^. The optimal 

initial state in this case is any pure state p^"^ with 
the reduced state equals to p, the simplest choice is 

| 00 ) 


v+ 




1+V+' ' V 1+V+' 

The maximal quantum Fisher information for spon¬ 
taneous emission without ancillary systems can also be 
calculated by imposing the state p used in M(p) to be 

where a, b are complex numbers 


pure: p = 
and la 


1 ^ 


|ar 

a*b 


\b\^ 


= 1. Eq. (IA32I) then becomes 


(si -I- 52 )^ =|a|'* + + 2 [pcosdx -f (1 - p)] |a6p.(A36) 


M(p) 


Pile +PP22e*'2' VI - pp2ie*^''+'2’) \ 

VI - ppi 2 e“*^^+'^^ (1 - p)p22F^ J ' 


Taking the derivative with respect to |ap gives the ex¬ 
tremal point |ap = ^ corresponding to the minimum 
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value 

(si + S 2 )^ = 1 + ^{cosdx — 1 ) = 1 — -rjdx^ + 0{dx*). 

(We know this is the minimum value by comparing to 
|ap = 0.) This gives the Bures angle 

cos[B{K^, K^+dx)] = 1 - ^vdx'^ + 0{dx‘^). (A37) 

O 


Finally, we get the maximal quantum Fisher information 
without using ancilla to be lim^^-i-o g i-cos[-B(J^x,-R'x+dx)] _ 
rj. The optimal input state is any pure state with pn = 
P 22 = 5 , which can be achieved by fQj- 
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